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We consider the astrophysical reaction rates for radiative neutron capture reactions (n, 7) in the 
crust of a neutron star. The presence of degenerate neutrons at high densities (mainly in the inner 
crust) can drastically affect the reaction rates. Standard rates assuming a Maxwell-Boltzmann 
distribution for neutrons can underestimate the rates by several orders of magnitude. We derive 
simple analytical expressions for reaction rates at a variety of conditions with account for neutron 
degeneracy. We also discuss the plasma effects on the outgoing radiative transition channel in 
neutron radiative capture reactions and show that these effects can also increase the reaction rates 
by a few orders of magnitude. In addition, using detailed balance, we analyze the effects of neutron 
degeneracy and plasma physics on reverse (7, n) photodisintegration. We discuss the dependence of 
the reaction rates on temperature and neutron chemical potential and outline the efficiency of these 
reactions in the neutron star crust. 

PACS numbers: 25.60.Tv, 26.60.Gj 



I. INTRODUCTION 

Nuclear reactions in the atmosphere and the crust of 
accreting neutron stars affect important observational 
manifestations such as X-ray bursts and supcrbursts 
(e.g., Rcfs. [H-ll]) as well as deep crustal heating of neu- 
tron stars in X-ray transients (e.g., Rcfs. 0-Q])- In the 
vicinity of the neutron drip density (p~4x 10 11 g cm' 3 
for the cold-catalyzed crust and p~6x 10 11 g cm~ 3 for 
the accreted crust Q) and beyond in the inner crust the 
dense matter contains an increasing amount of free de- 
generate neutrons (see, e.g., Ref. [8[). Neutron capture 
and reverse reactions are important components of nu- 
clear burning under these conditions @. Standard ther- 
monuclear neutron capture rates, which are used in re- 
action network simulations of nucleosynthesis in stars or 
supernova explosions, are obtained (e.g., Ref. [lOj]), as- 
suming the classical Maxwell-Boltzmann distribution of 
neutrons. However, free neutrons in the neutron star 
crust can be degenerate, in particular when the density 
exceeds the neutron drip point 0. For instance, ground- 
state (cold-catalyzed) matter at p — 6.2 x 10 12 g cm -3 has 
a neutron Fermi energy of « 2.6 MeV Consequently 
neutron degeneracy needs to be taken into account for 
neutron capture rates under such conditions. 

In addition, the dense stellar plasma of the neutron 
star crust strongly affects emission, absorption, and prop- 
agation of photons fl2j and therefore modifies radiative 
capture and photodisintegration reactions, like (n, 7) and 
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(7, n). Because of the high density, the electron plasma 
frequency uj p can be of the order of or higher than char- 
acteristic frequencies of radiative transitions in nuclei. 
Under these conditions, well-defined elementary electro- 
magnetic excitations (photons or plasmons) become ei- 
ther suppressed or forbidden (e.g., Ref. [13|) although 
radiative transitions are not suppressed because they can 
be realized by emission (or absorption) of excess energy 
to (from) the plasma as a collective system [l2|. These 
plasma physics effects can be important since they may 
enhance the radiative transition strength. 

In Sec. |n] we discuss the effects of neutron degener- 
acy on (n, 7) radiative neutron capture reactions in dense 
matter. In Sec. IHII we analyze plasma effects on the out- 
going radiative transition channel of (n, 7) reactions. In 
Sec. IIVI we consider the same neutron degeneracy and 
plasma physics effects on inverse (7, n) photodisintegra- 
tion reactions. We discuss our results in Sec. M and sum- 
marize them in Sec. I VII For brevity, we use the units in 
which the Boltzmann constant Hb = 1. 



II. REACTION RATES FOR DEGENERATE 
NEUTRONS 

We start with the outline of the (n, 7) radiative capture 
rates in stellar environments (e.g., Ref. [HI). Let the 
cross section cr a b(E) refer to the reaction X^ + n — > 
y( b ) + j t where a and b label different energy levels of a 
target nucleus X and a resultant nucleus Y , respectively, 
and E is the center-of-mass energy of the reactants. In 
stellar matter at local thermodynamic equilibrium, the 
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total cross section <r*(E) of the reaction X + n —> Y + 
7 includes neutron capture on the ground state and all 
thermally populated states, 



a*(E) 



Ea9aexp(-E^/T)j: b Va b (E) 
J2 a g a exp(-E¥/T 



(1) 



where E^ is the energy of level a and g a is its statistical 
weight. The summation in the denominator normalizes 
the distribution of target nuclei over the energy levels a 
(gives the internal partition function of the target nu- 
cleus). Asterisk * means that thermally excited nuclear 
levels in stellar matter arc included. 

The astrophysical reaction rate contains the average 
(cr*v) of the total cross section with the energy distribu- 
tion f{E) of the interacting particles. For nonrelativistic 
reactants (considered in this paper) the collision energy 
is E = /j.v 2 /2, where /i is the reduced mass (very close 
to the neutron mass m n ), and v is the relative velocity 
of a neutron with respect to nucleus at large separations. 
Then 



^ v) = {7^[ Ea * {E)mAE ' 

where Af is the normalization factor 

Af = VE f(E) dE. 



(2) 



(3) 



We call (a*v) [cm 3 s _1 ] the reaction rate coefficient. The 
rate itself (for instance, per unit volume, cm~ 3 s -1 ) is 
nxn n (a*v) , where n n and nx are number densities of 
neutrons and reacting nuclei, respectively. 

Astrophysical reaction rates at typical stellar temper- 
atures are based on a Maxwcll-Boltzmann distribution of 
the particles, Jmb(E) = exp(—E/T). At high densities in 
the neutron star crust, neutrons can become degenerate, 
which modifies the reaction rate. At these conditions the 
nuclei are not freely moving particles but are confined in 
a strongly coupled Coulomb liquid or a Coulomb crystal 
(e.g., Ref. @). Because the neutrons arc much lighter 
than the nuclei, the energy distribution function f(E) in 
Eq. ([2]) can be approximated by a Fermi-Dirac distribu- 
tion /fd(-E) = [1 + cxp((E — n n )/T)]~ , where fi n is the 
neutron chemical potential. In this approximation we ne- 
glect recoil effects and nucleus motion. In the inner crust 
of the neutron star the nuclei can be bulky and occupy 
a non-negligible fraction of volume fl4| . Here we employ 
the model of a free neutron gas with local number density 
n n which occupies the space between the nuclei. Though 
this model is rather accurate near the neutron drip point , 
it becomes less accurate at higher densities where free 
neutrons constitute a strongly interacting Fermi liquid 

Let the average (a* w)mb be obtained with the 
Maxwcll-Boltzmann distribution and (cr* u)fd be calcu- 
lated with the Fermi-Dirac distribution. Many calculated 



reaction rate coefficients (i7*t;)mb for neutron capture re- 
actions are available in the literature (e.g., Ref. [l5[ and 
references therein). In a neutron star crust (<7*i>)fd de- 
pends on T and \x n (or, equivalently, on T and n n ). For 
practical applications we introduce the ratio 



Rr, 



\Q V)FD 
(cr*v) MB 



(4) 



These ratios are easier to calculate and approximate than 
(cr*t;)FD; the derivation of these ratios ise the main sub- 
ject of the present paper. 

Generally, accurate calculations of (c*i;)fd and R n 
require the cross sections a a b{E) obtained from exper- 
iment or nuclear reaction codes. Detailed calculations of 
(<7*d)fd would be a valuable project for the future. Here 
we restrict ourselves to a simplified approach. It allows 
us to demonstrate the importance of the effects of neu- 
tron degeneracy and is sufficiently accurate for a wide 
range of temperatures and densities. First, we neglect 
the contribution of thermally excited states (setting thus 
a*(E) = cr(E)). This is a valid approach if the energy of 
the first excited level of the target nucleus is higher than 
the temperature in the neutron star crust (T < 2 x 10 9 
K« 0.2 MeV). For threshold (endothermic) reactions, it 
should also be higher than the reaction threshold Eq. 
Second, we note that the reaction rates at low tempera- 
tures correspond to cross sections o{E), which are char- 
acterized by typical power-law behavior that we therefore 
adopt in our analysis: 



a(E) = cr a (E- EoY at E < E < E n 



(5) 



Here, v is a power-law index, a a is a normalization con- 
stant, Eq is a reaction threshold (with Eq = for exother- 
mic reactions), and E max is the maximum energy up to 
which the approximation ([5]) holds. We treat Eq, a a , v, 
and E max as input parameters. For a given reaction, they 
can be adopted from a nuclear database or calculated us- 
ing a nuclear reaction code. In our approximation, the 
factor R n depends on T, fi n , v, and Eq; the parameter 
<t q cancels out in the ratio E max is required to check 
the validity of the calculated R n for given conditions. 

The reaction rates are strongly affected by the neu- 
tron energy distribution. Typical energies of nondegen- 
erate neutrons are E < T. In a strongly degenerate 
gas (/i„ 3> T) the majority of neutrons belong to the 
Fermi sea and have much higher energies T <C E < /j, n . 
In this case, there is also a smaller (but non-negligible) 
amount of neutrons, with energies above fi n : fi n < E < 
fjt n + T. Their distribution /fd(-E) ~ exp((^„ — E)/T) = 
exp(^„/T) /mb(E) is close to Maxwellian and represents 
the Maxwellian tail of the Fermi-Dirac distribution. In 
the following we demonstrate that these different energy 
ranges of /fd (E) correspond to different neutron capture 
regimes. 

For all nuclear reactions shown in the paper we use 
the cross sections obtained with the statistical model 
Hauser-Feshbach (HF) code TALYS-1.2 [Hj]. Statistical 
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model theory [17[ uses the concept of averaged trans- 
mission coefficients to describe the formation and sub- 
sequent decay of a compound nucleus formed after a 
projectile impinges on a target nucleus. In this sce- 
nario the reaction sequence for neutron capture becomes 
+ n — > C* — > + 7, where C* is a compound nu- 
cleus with many closely spaced energy levels (high level 
density). For the neutron star crust conditions the in- 
cident neutron has rather low energy and the primary 
reaction mechanism is dominated by compound forma- 
tion. The partial cross section a a b (E) in the HF model is 
written as the sum over levels c (specified by energy E Cl 
spin J, and parity 7r) of the compound nucleus 



each allowed channel is obtained by solving the 
Schrddinger equation with an optical potential for the 
neutron-nucleus interaction. The 7-transmission coeffi- 
cient is calculated for a giant dipole resonance (El+Ml) 
approximated b y a sing le Lorentzian or by a combination 
of Lorentizains [18l - l20l ] . The sum of these contributions 
determines the 7-ray strength function. Both the neu- 
tron and 7-ray transmission coefficients must be calcu- 
lated for all accessible states. In practice there is a huge 
number of levels, the vast majority of which are exper- 
imentally unexplored. For very neutron-rich nuclei near 
the drip line the level density may be much smaller and 
the applicability of the HF model may be questionable. 



9c 



' n,a 7,6 



9n9a T t C c 



(0) 



In this case, k is the wave number of an incident neutron, 
T„ a and T^ b are partial transmission coefficients, and 
7^Q t is the total transmission coefficient of the compound 
nucleus in a level c. The latter quantity, T^ ot = b 7~ob' 
gives the total width of the c level as the sum over all 
available outgoing reaction channels o — n, 7, etc. and 
over levels b of the final nucleus. Note, that the sum 
includes compound elastic scattering (when final states 
are the same as the initial ones). 

The individual neutron transmission coefficient for 



Other necessary ingredients for a HF calculation in- 
clude the choice of level density, optical model, 7-ray 
strength function, and mass model to predict the reac- 
tion Q value. The reaction cross sections presented here 
do not include pre-equilibrium effects. The cross sections 
are calculated on the basis of Q values derived from the 
Hartree-Fock-Bogoliubov mass model HFB-17 [2l|. The 



level densities are obtained from the microscopic model of 
Ref. [22|, and the El 7-ray strength function is based on 
quasiparticle random-phasc-approximation calculations, 
folded with a simple Lorentzian [23|. The neutron opti- 
cal potential is supplied by the global parametrizations 
of Ref. [H. 
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FIG. 1. Cross sections a of neutron capture on 39 Mg (panel 
(a)) and 46 Mg (panel (b)) plotted (left vertical axis) in double 
logarithmic scale as a function of (E — Eq)/T at T = 0.1 MeV. 
In both panels we show also (right vertical axis) the filter 
functions Ef(E)/J\f for the three cases: Maxwell-Boltzmann 
distribution (dashed lines) and Fermi-Dirac distribution with 
fi„ — 1 MeV (dash-dotted lines) and /u n = 5 MeV (dotted 
lines). See text for details. 



Figure [T] shows the reaction cross sections (left vertical scales) for two neutron capture reactions, on 39 Mg (panel 
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(a)) and 46 Mg (panel (b)). The 39 Mg(n, 7 ) 40 Mg reac- 
tion is exothermic (E = 0), while the 46 Mg(n, 7) 47 Mg 
is endothermic (Eq — 4.06 McV). For a better visual- 
ization of the approximation (JS|), u{E) is shown as a 
function of (E — Eq)/T on a double logarithmic scale; 
the temperature is taken to be T = 0.1 MeV. The lin- 
ear segments of the curves clearly indicate the power-law 
behavior of cr(E) at low E. Power-law indices and max- 
imum energies are v = —0.6 and E max w 0.1 MeV for 
neutron capture on 39 Mg; v = 3.5 and E max ps 6 MeV 
for neutron capture on 46 Mg. This power-law behav- 
ior at low energies is typical for (n, 7) reactions. Fig- 
ure [T] also shows (right vertical scales) the so-called filter 
functions Ef(E)/Af, which enter the integrand of (j2|) 
along with a*(E). The dashed lines in both panels cor- 
respond to the Maxwell-Boltzmann distributions of neu- 
trons. We see that the power-law approximation is suf- 
ficient for calculating (<tv)mb m both cases (at T = 0.1 
MeV). The dash-dotted and dotted lines in Fig. [1] rep- 
resent the Fermi-Dirac distribution with fi n = 1 MeV 
and 5 MeV, respectively. We see that for /i„ = 1 MeV 
and the 46 Mg target the power-law approximation is def- 
initely valid, while for fi n = 5 McV it is less accurate. 
For neutron capture on the 39 Mg nucleus, the power-law 
approximation is inaccurate at both values of fi n . 

It is easy to see that the power-law approximation is 
valid as long as max(£'o, fi n )+T < E max . In this approxi- 
mation, the factor R n in Eq. ((4]) is calculated analytically. 
By introducing dimensionless parameters y = \i n jT and 
xq = Eq/T, we obtain, 



Rn 



expx {v + l)F v+ i(y - x ) + x F u {y - x Q ) 



x + u+l Fi/2(y) 
where F v (y) is a Fermi-Dirac integral 



My) 



1 



x v dx 



T{v+\)J l + exp(x-y) 



(7) 



(8) 



and T(y + 1) is the Eulcr gamma-function. 

Equation ([7]) is the main result of our consideration. 
We expect that this factor is sufficient to correct the 
reaction rate for neutron degeneracy in many cases of 
practical importance. Let us analyze the limiting cases. 
For this purpose we use asymptotes of the Fermi-Dirac 
integrals given in the Appendix. 

First consider a threshold reaction with typical neutron 
energies well below the threshold. In this case xq — y 1 
(E - fJ, n > T), and Eq. (JTJ) becomes 



Rn — 



cxp y 



(9) 



If, in addition, neutrons are strongly degenerate ((/> 1), 
Eq. ([9]) is further simplified, 



Rn 



-^—y 3/2 expy. 



The factor R n in Eqs. (j9)) and (p~0|) becomes a function 
of only one parameter, y = H n /T\ it is independent of 
v. This indicates that Eqs. dU and (fTU|) arc valid for 
any cross sections cr(E) of threshold reactions, not only 
for ones with power-law behavior. Indeed, in the limit 
ill - 1/ > 1 the majority of neutrons have energies E < 
Eq; these neutrons cannot overcome the reaction thresh- 
old and cannot be captured by nuclei. The reaction pro- 
ceeds owing to a small amount of high-energy (suprathcr- 
mal) neutrons with E > Eq. Recall that for any neutron 
degeneracy, their distribution function in Eq. ([5]) is ac- 
tually Maxwellian, /fd(-E) w exp (/j, n /T) /mb(-E)- Then 
the cross sections in the nominator and denominator of 
Eq. (jlj are integrated with the same function /mb(E); 
equal integrals for any a(E) cancel out and do not affect 
the ratio R n . 

The ratio R n in Eq. fjlO|) shows a sharp exponen- 
tial y-dependence at strong neutron degeneracy. This 
means that neutron degeneracy exponentially enhances 
the rates of threshold reactions (by increasing the amount 
of high-energy neutrons). 

Another limiting case is the one of an (endo- or 
exothermic) reaction with strongly degenerate neutrons 
whose typical energies are above the reaction threshold. 
In this limit y — xq 3> 1 (/x n — Eq 3> T), the reaction 
is driven by numerous energetic Fermi-sea neutrons and 
becomes fast, with 



Rn 



3^7? CXp .To (y - XqY Xq + (y + \)'IJ 

4T(j/ + 3) y 3 / 2 X0 + 1/ + I 



(11) 



In particular, for a threshold reaction with strongly 
degenerate neutrons at Eq 3> T and /x„ 3> T, R n is 
given by Eq. flTDJl for E - /j, n > T and by Eq. (JTTJ) for 
/i„ — Eq ^> T. These two asymptotes nearly match each 
other at \fx n — Eq | ~ T providing an accurate description 
of R n in a wide range of /i n except for the narrow interval 
iMn - E \ < T that should be described by Eq. 0. 

At /j, n > Eq the dependence of R n on the neutron 
degeneracy parameter y is much weaker than at /i„ < 
Eq. In the range of T < n n - Eq < Eq, Eq. (JTTJ) can 
be simplified by setting y = xq everywhere but in (y — 



(10) 



xq) u+ ; this gives R n oc (//„ — E ) u+1 . For higher fj, n ^> 
Eq we have R n oc [i n +1 ^ 2 ■ 

In order to calculate R n in intermediate cases from 
Eq. 0, accurate expressions for Fermi-Dirac integrals 
arc necessary. These integrals have been extensively 
studied in the literature, especially in the field of semi- 
conductor physics and astrophysics (e.g., Ref. [lEHl} and 
references therein). There are several very accurate ap- 
proximations for particular integer and half-integer val- 
ues of v. An analytic approximation that can be used 
for any v in the range —l<v<A was constructed 
by Aymerich-Humet et al. |27| . It accurately repro- 
duces the limits of y — >• ±00, and its relative error at 
— 1/2 < v < 5/2 does not exceed 1.2%. For convenience, 
we present this approximation in the Appendix. 
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III. PLASMA EFFECTS 

In addition to the effects of neutron degeneracy, the 
rates of (n, 7) reactions in dense matter are influenced by 
electron plasma effects. Under typical conditions in the 
neutron star crust, the electrons behave as weakly inter- 
acting, strongly degenerate and ultra-relativistic parti- 
cles [8[ . The importance of plasma effects is characterized 
by the electron plasma frequency oj p = ^/47re 2 n e /m*, 
where e is the electron charge, n e is the electron number 
density, and m* = \i e j<? is the effective electron mass at 
the Fermi surface (fi e being the electron chemical poten- 
tial). The plasma effects modify the radiative transition 
in the exit channel. The plasma effects are strong when 
the frequency u of a radiative transition becomes com- 
parable to u>p. In particular, no well-defined electromag- 
netic excitations (photons or plasmons) can propagate 
at ui < uj p (e.g., Ref. [IH). In this case the radiative 
transition cannot occur through the emission of a real 
photon or plasmon. However, it can occur through a 
direct transfer of the excess energy ftiu to plasma elec- 
trons via collision-free collective electromagnetic interac- 
tions. For a degenerate electron plasma in the neutron 
star crust, this effect was considered in Ref. It does 
not suppress but rather enhances radiative transitions 
at uj < Lu p . The enhancement factor 1Z 7 (lu) depends on 
the radiative transition type (electric or magnetic), on 
the transition multiplicity £ — 1,2,..., and on the lu/uj p 
ratio. In the low-frequency limit, 



behaves as 7?. 7 (w) ~ (uip/u)) 21 , while at w » Wj, one has 
lZ 1 {w) — > 1, meaning that the plasma effects become less 
important and the standard regime of emission of real 
photons is restored. 

The plasma effects modify the reaction cross section 
<j(E) and the reaction rate coefficient (cry). In analogy 
with Eq. (0J we can formally introduce the total correc- 
tion factor 



w « w pi this factor 



R 



72, pi 



(< J P \ V ) 



FD 



ffW MB 



(12) 



which takes into account both neutron degeneracy and 
plasma effects; cr{E) and & P i(E) are the cross sections 
excluding and including the plasma effects, respectively. 
It is convenient to write 



(13) 



where R n takes into account neutron degeneracy alone as 
discussed in Sec. HIl and R p i accounts for plasma effects 
(in the presence of neutron degeneracy). 

The inclusion of plasma effects in (n,, 7) reaction rates 
is not straightforward. As discussed in Sec. HH the rates 
of interest are usually calculated [l(| in the framework 
of the HF statistical model [TtJ ■ The radiative transition 
coefficients need to be modified in the numerator and 
denominator of Eq. ([5]) . The problem is further compli- 
cated by summing over thermally excited nuclear levels 
b in Eq. ([T]). For different levels b the energy and type 



of radiative transition can be different (implying differ- 
ent plasma modifications). A correct inclusion of the 
plasma effects is therefore a complicated computational 
project, which is beyond the scope of our paper. Here, 
we present a simplified approach which demonstrates the 
importance of the plasma effects. It is based on the as- 
sumption that the radiative transmission coefficients Tj 
are much smaller than other contributions to the total 
transmission coefficient 7tot in Eq. (|6|). In that case the 
total transmission coefficient 7tot is independent of the 
plasma effects. The second assumption is that the radia- 
tive exit channel is represented by a single El radiative 
transition to the ground state, so that no summation over 
excited states b is required. 

These requirements are realized for threshold (n, 7) re- 
actions involving degenerate neutrons where fi n is below 
or slightly above Eg. In this case the typical energy of 
interacting neutrons is not much higher than Eq, which 
reduces the radiative decay of the compound states to a 
single (not multiple) low-energy radiative transition. 

Under these assumptions we can describe the plasma 
effects by multiplying the neutron capture cross section 
<j(E) by the factor 1Z ie \(uj), with huj = E — Eq, which 
describes the enhancement of the radiative transition by 
the plasma effects. The latter factor was calculated in 
Ref. [lH and fitted by a simple analytical expression 
(Eq. (33) in HU) with an error of about 1%. The mod- 
ified cross section must be integrated over E to obtain 
the reaction rate. To simplify the integration we suggest 
using the approximation 



TZ jel {uj) w 1 + 3.03^)' 



(14) 



(jj p , with the maximum 



Although it is rather crude at uj 
error 17% at to — uj Pi it reproduces the correct asymptotic 
behavior for u> 3> ui p and uj <^ lu p . The deviation does 
not exceed 5% outside the region of 0.6 < uj/uj p < 1.03. 
The advantage of using Eq. (fl"4")l is that it allows an ana- 
lytic integration of the correction to the reaction rate in 
the same power-law approximation for the cross section 
as used in Sec. Ill Because 1Z ie i{uj) is integrated, the ap- 
proximation errors cancel out, leaving us with a rather 
accurate result. The correction factor (|T2"j) to the reac- 
tion rate including both neutron degeneracy and electron 
plasma effects now becomes 



3.03 



v(y -\)\T 



x 



v -1 



x 



1 



Rir 2) + rL i,) - 



(15) 

Here, R^f 1 is the factor given by Eq. ([7]) for a power-law 
index v. 

In the limit Eq — [i n ^> T for a threshold reaction 
with Eq ^ T, the expression (fl~5j) is further simplified. 
Introducing the correction factor R p \ due to the plasma 
effects in accordance with Eq. (|13p. we have 



Rpi = 1 + 



3.03 



FlJUJj 

~~T 



(16) 
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In this limit, in addition to a strong exponential enhance- 
ment of the reaction rate due to neutron degeneracy (Sec. 
Hi]) , there is a smaller but significant plasma enhance- 
ment. In the opposite limit of fj, n — Eq 3> T (but still for 
a single radiative transition) we obtain 



Rpi = 1 + 3.03 



(i/ + 2)(i/+ 1) fUJ p \ 2 E + (u- 1>„ 



!/(!/-!) 



+ (v+lK 



(17) 

This factor is temperature independent because the typ- 
ical transition energy is now nearly fixed by /z n and Eq, 

tlLJ RJ H n - Eq. 



IV. RATES OF INVERSE REACTIONS 

If the rate of a forward X(n,j)Y reaction is known, 
then the rate for an inverse reaction Y(j, n)X can be 
determined from the detailed balance principle, 



where ; n y°^ > ano - n " q; are number densities of nu- 
clei X, Y , and neutrons, respectively, in statistical equi- 
librium between forward and inverse reactions. A 7 [s _1 ] 
specifies the rate of the inverse (photodisintegration) re- 
action (which is ny\ 7 , cm -3 s _1 ). Usually, this reaction 
involves only photons 7 but in our case it also can involve 
more complicated excitations associated with the elec- 
tromagnetic field and plasma electrons (Sec. IIII[) , which 
are assumed to be in thermal equilibrium; their effective 
number density is included in A 7 . The equilibrium num- 
ber densities of the nuclei should satisfy the condition of 
chemical equilibrium 



>q) 



(19) 



where fix and \xy are the chemical potentials of the nu- 
clei X and Y, respectively. Traditionally one assumes 
ideal nondegenerate gas conditions for the nuclei and neu- 
trons in order to relate their equilibrium number densities 
and chemical potentials. In this approximation Eq. (TH 
yields the well-known relation 



A x m n T 
2nh 2 A Y 



3/2 



cxp ^-^ ) 



where Ax and Ay = Ax + 1 are mass numbers of nu- 
clei X and Y, respectively, while Zx and Zy are their 
individual internal partition functions 



E 



9a exp 



T 



(21) 



The corresponding function for the neutrons is Z n = 2. 

Now we should modify Eq. ([277)) to account for neutron 
degeneracy and strong Coulomb coupling of the nuclei. 



Strong Coulomb coupling prevents treating the plasma 
of atomic nuclei as an ideal gas (it becomes Coulomb liq- 
uid or crystal; e.g., Ref. ||). A strongly coupled multi- 
component plasma of charged particles satisfies (to a very 
high accuracy) the linear mixing rule according to which 
the main Coulomb thermodynamic quantities (like mean 
Coulomb energy, etc.) can be presented as sums of quan- 
tities for individual ions. Coulomb coupling of an individ- 
ual atomic nucleus X = (Ax, Zx) in this plasma is de- 
scribed by the parameter Tx = Z x e 2 /(axT), where ax 

1/3 

is the ion sphere radius defined as ax = [3Zx/(4-7m e )] 
This allows one to treat a strongly coupled system of 
atomic nuclei as an ensemble of weakly interacting ion 
spheres. This approximation is well known in the physics 
of strongly coupled Coulomb plasmas Q. In this case 
Eq. (|20[) remains the same but the internal partition func- 
tion for each nucleus has to be multiplied by its individ- 

(c) (o 
ual Coulomb partition function Z x . A function Z x 

depends only on one parameter, Tx, which, in turn, is 

determined by the nuclear charge number Zx- Because 



(18) in our case Zx = Zy , we have Z 



(C) 
X 



Z^, and the 

Coulomb corrections for the nuclei X and Y compensate 
for each other in Eq. (f2T)j) . 

Neutron degeneracy can be included in Eq. (f2"Uj) by im- 
plying the correct relation between the neutron number 
density and its chemical potential. It is easy to show that 
for this purpose it is sufficient to multiply the right-hand 
side of Eq. (J20j) by 



i? rvs = exp(-y)J r 1/2 (y). 



(22) 



The ratio of the photodisintegration rates for the Fcrmi- 
Dirac and the Maxwell-Boltzmann distributions of neu- 
trons then becomes 



(A 7 )fd 
(A 7 )mb 



R\ = R IVS R„ 



(23) 



If plasma effects are included, then R n must be replaced 
by Rn, p v 

Various asymptotes for R\ are readily obtained from 
those for R n , Eqs. (|9|)- (fTT|) . In particular, for the case 
of a threshold neutron capture reaction with the neutron 
chemical potential well under the threshold, Eo—jj, n ^> T, 
we obtain (neglecting plasma effects) R\ = 1 for any de- 
pendence of a on E. The inverse reaction is not affected 
by neutron degeneracy which is quite natural. The ef- 
fect of neutron degeneracy on an (7, n) reaction consists 
of Pauli blocking of the emitted neutrons. However, in 
our case these neutrons have energies E = Q + huj 3> 
above the Fermi level, where the blocking does not occur. 

In the opposite limit, fi„ — Eq 3> T, the neutrons, 
emitted in the reverse reaction, have low energies and 
are strongly blocked by the Fermi sea neutrons; this ex- 
ponentially suppresses the inverse reaction rate: 

(y - xor +1 [xp + (y + l)y] 

R\ = — 7 ; TTTTv — 7T\ — ex P - y). 24) 

(x + v + 1) Y[y + 3) 
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Note that recently Mathews et al. [28| suggested mod- 
ifying the detailed balance equation (f2T))) by taking into 
account the quantum corrections due to induced pho- 
ton effects in the photodisintegration rate coefficient A 7 . 
They point out that while calculating A 7 one usually em- 
ploys the Maxwellian distribution of photons instead of 
the Planck distribution fp\ — (cxp(£' 7 /T) — Using 
the Planck distribution, they corrected A 7 and concluded 
that one should also correct the detailed-balance ratio 
A 7 /(ot>)mb- However, in this latter conclusion the au- 
thors erroneously neglected the same corrections in the 
rate coefficient (ctw)mb of the forward reaction. Specifi- 
cally, they did not include an extra factor (f + /pi(i? 7 )) 
(with E 1 = E + Q) under the integral in their Eq. (6) 
(similar to Eq. ([2]) in the present paper) to account for 
the induced emission. If that factor would have been in- 
troduced, the detailed balance ratio would be naturally 
unaffected by quantum corrections. In our analysis we 
neglect such corrections in both forward and reverse re- 
action rates, because they are generally small [28j . 



V. DISCUSSION 




1 1 r— I 

-2 -1 



log lo r[MeV] 

FIG. 2. Factor R n for the 46 Mg(n, 7) 47 Mg reaction as a 
function of temperature for fj, n = 1, 3, and 5 MeV. 

Let us illustrate the obtained results. First consider 
the effects of neutron degeneracy on threshold reactions 
neglecting plasma physics effects. Figure [2] shows R n 
as a function of T for the 46 Mg(n, 7) 47 Mg reaction at 
fx n = 1, 3, and 5 MeV. This reaction has a rather high 
threshold, E = 4.06 MeV. We see a strong increase of 
the reaction rate with growing fi n . For fi n — 1 and 3 MeV 
the factor R n is well described with the power-law ap- 
proximation by Eq. ([9J. For ^„ = 5 MeV the difference 
to the power-law approximation becomes noticeable, but 
not on the logarithmic scale of Fig. [2] 

To illustrate the possible deviations from the power- 
law approximation, in Fig. [3] we plot R n versus [i n at 
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FIG. 3. Factor R n for the 62 Ca(n,7) 63 Ca reaction as a func- 
tion of fi n at T = 0.1 MeV. Vertical dotted line indicates the 
reaction threshold. The inset shows the reaction cross section. 
Dashed lines in the figure and inset refer to the power-law ap- 
proximation. 



T = 0.1 MeV (the main figure) and a(E) (the inset) for 
the 62 Ca(n, 7 ) 63 Ca reaction (E = -Q = 0.78 MeV). 
The solid lines are obtained by numerical calculations, 
while the dashed lines are the results of the power-law 
approximation. It can be seen that the latter approxi- 
mation accurately describes a{E) up to £ ma s « 10 MeV. 
Accordingly, Eq. (JT]) closely reproduces the dependence 
of R n on fi, n at [i n < 10 MeV; the exponential asymptote 
pTJl) is valid at [i n < 1.5 MeV; the power- law asymp- 
tote ([11]) works well at 1.5 < fi n <10 MeV. At [i n > 10 
MeV the power-law approximation becomes inaccurate 
because at such fi n the reaction rate is affected by the 
high-energy segment of cr(E) where the power-law is in- 
valid. 

The situation is different with exothermic reac- 
tions. Figure 0] shows the calculated factor R n for the 
39 Mg(n, 7) 40 Mg reaction. The reaction is exothermic 
(E =0, Q = 1.4 MeV); its cross section is plotted in the 
left panel of Fig. [TJ The solid, dashed, and dot-dashed 
lines in Fig. |4] are calculated for /i„ = 0.2, 1, and 5 MeV, 
respectively. Thick lines are obtained by integration of 
numerically calculated cross sections; thin lines are ob- 
tained using the power-law approximation. We see that 
the effect of neutron degeneracy is much weaker than for 
threshold reactions; the factor i? n stays ~ 1. Notice that 
for /j, n = 0.2 and 1 MeV the Fermi-Dirac averaged rate 
is smaller than the Maxwell-Boltzmann rate (R n < 1). 
The power-law approximation breaks down even for small 
values of /i„ {o{E) starts to deviate from power-law at 
E < 0.2 MeV, see Fig.©. 

Now let us discuss the impact of plasma effects on neu- 
tron capture reactions starting with the threshold reac- 
tions. In all cases the plasma effects cause an additional 
enhancement of the reaction rates. Figure [5] gives R p \ for 
the 38 Mg(n, 7) 39 Mg reaction as a function of fi n at the 
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T [MeV] 

FIG. 4. Factor R n for the 39 Mg(n, 7) 40 Mg reaction as a 
function of temperature at \i n — 0.2, 1, and 5 MeV (solid, 
dashed, and dot-dashed lines, respectively). Shown are re- 
sults of numerical integration of calculated cross section a(E) 
(thick lines) and results of the power-law approximation (thin 
lines). 
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FIG. 5. Factor R p i vs fx„ for the 38 Mg(n, 7) 39 Mg reaction in 
matter with the electron plasma frequency ftoj p = 1.5 MeV. 
Solid, dashed, and dash-dotted curves are plotted for T = 0.1, 
0.5, and 1 MeV, respectively. Thick curves are obtained using 
computed cj{E) and accurate approximation of TZ le i(uj); thin 
curves represent the approximation (|15|l . Vertical dotted and 
dash-dotted lines indicate the reaction threshold and plasma 
frequency, respectively. 



the typical energy released in the radiative transition is 
ftio ~ T, while at \i n > Eq it is fku m \i n — Eq. Therefore, 
at /i„ < Eq the factor R p \ is mostly independent of /z n , 
while at /j,„ > Eq it is independent of T (cf. Eqs. p^|) 
and (|17p ). A significant plasma enhancement can be 
observed at fi n < Eq. It is most visible for the solid 
curves (in those segments where T and /i„ are smaller 
than ojp). The plasma enhancement is expected to be 
especially pronounced in those reactions which occur in 
dense plasma and are accompanied by small radiative 
energy release. 

The difference between the thick (more accurate) and 
thin (less accurate) lines in Fig. [5] is small for T = 
0.1 MeV but is higher for T = 0.5 and 1 MeV. This 
is solely due to the breakdown of the power-law ap- 
proximation of <j{E) for the 38 Mg(n, 7) 39 Mg reaction at 
E > 1.5 MeV. We have checked that using Eq. (H]) in- 
stead of the more accurate approximation to TZ-yei (w) 
from Rcf. [12] does not significantly change the results. 

For exothermic reactions the simple model used in 
Sec. IIIII is generally inapplicable. We expect that plasma 
effects on these reactions give R p i ~ 1. 
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FIG. 6. Lines of constant effective times r = 10 9 , 1, 10 9 s 
of 46 Mg burning in the 46 Mg(n, 7) 47 Mg reaction on the T-fi n 
plane. We present calculations using the Maxwell-Boltzmann 
neutron capture rate (dashed lines), the rate corrected for 
neutron degeneracy (solid lines), and the rate corrected for 
degeneracy and plasma physics effects (dash-dotted lines). 
Vertical dotted line positions the reaction threshold Eq = 
4.06 MeV. The double-dot-dashed line corresponds to equal 
amounts of 46 Mg and 47 Mg nuclei, assuming statistical equi- 
librium with respect to neutron capture and emission reac- 
tions. 



electron plasma frequency hio p = 1.5 MeV. The curves 
are plotted for T = 0.1 (solid lines), 0.5 (dashed lines), 
and 1 MeV (dash-dotted lines). Thick curves are cal- 
culated with numerically determined cross section and 
an accurate plasma enhancement factor 7£ 7e i(uj) from 
Ref. [13 • Thin curves are given by simplified Eq. (fT31) . 
All curves change their behavior when fi n reaches the re- 
action threshold Eq = —Q = 0.5 MeV. For fi n < Eq. 



Finally, in order to illustrate the efficiency of neu- 
tron captures, in Fig. [6] we plot the lines of constant 
characteristic burning times t of 46 Mg nuclei in the 
46 Mg(n,7) 47 Mg reaction. The plot is made in the T- 
fi n plane for the three values of r = 10 -9 , 1, 10 9 s. The 
characteristic burning time is defined as 

T- 1 =n n (<rv). (25) 
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In the region above and to the right of each line burning 
is faster than on the line (t is smaller); in the region 
below and to the left burning is slower. Lines of different 
types are calculated using different reaction rates. The 
dashed lines are for non-degenerate neutrons, the solid 
lines take into account neutron degeneracy, and the dot- 
dashed lines take into account neutron degeneracy and 
plasma effects (assuming fku p = 1.5 MeV). Neutrons are 
essentially non-degenerate in the region of fx n < (in 
the left part of the fi n — T plane). They are strongly 
degenerate in the region of /i„ > T. 

As long as neutrons are non-degenerate, the solid and 
dashed lines naturally coincide. When neutron degen- 
eracy sets in, it intensifies neutron captures. Then the 
solid lines (in sharp contrast with the dashed ones) bend 
and become nearly vertical. In this regime neutron de- 
generacy is vitally important. The most remarkable ef- 
fect occurs in the vicinity of the threshold (fj,„ = Eq = 
4.06 MeV, shown by the dotted vertical line). We see 
that the solid lines drop off to zero temperature almost 
immediately after p n exceeds the reaction threshold. In 
this case the reaction is driven by the Fermi sea of de- 
generate neutrons and becomes extremely fast. 

The purpose of Fig. [5] is primarily to illustrate the ef- 
ficiency of the 46 Mg— s- 47 Mg transformation with respect 
to neutron capture, neglecting other reactions (reverse 
reaction, beta captures, and fusion reactions). Its main 
and natural result is that 46 Mg cannot survive for a long 
time against neutron capture when \i n exceeds Eq. 

According to Fig. [51 the 46 Mg— > 47 Mg reaction actu- 
ally occurs in a narrow strip on the T — p n plane that 
is confined between the lines r = 10~ 9 s and r = 10 9 
s. Below and to the left of the r = 10 9 line, r is very 
large; there will be plenty of 46 Mg nuclei which are very 
inefficient neutron absorbers. Above and to the right of 
the r = 10 -9 s line, r is extremely short; all 46 Mg nuclei 
are transformed into 47 Mg. The shape of this "burning" 
strip is similar to that for fusion reactions; see, for in- 
stance, Fig. 4 of Ref. (2{|, which gives the strip for the 
12 C+ 12 C fusion reaction in the T — p plane. The bend of 
the carbon burning r=const lines at high densities p is 
produced due to the transition from thermonuclear car- 
bon fusion to pycnonuclear one. It greatly resembles the 
bend of the r=const curves with the growth of /i n due to 
the effects of neutron degeneracy in our Fig. [5] 

In addition, the double-dot-dashed curve in Fig. [J5] is 
the line representing equal amounts of 46 Mg and 47 Mg 
nuclei (n 1 ^ = riy q ^), assuming statistical equilibrium 
with respect to the neutron capture and reverse reactions. 
In statistical equilibrium, the matter (in our simplified 
model) would mainly contain 46 Mg nuclei below this line 
and 47 Mg nuclei above this line. 



VI. CONCLUSIONS 

We have considered neutron captures (n, 7) in dense 
stellar matter, taking into account the effects of neutron 



degeneracy and plasma physics. The effects of neutron 
degeneracy increase the amount of high-energy neutrons 
and mainly enhance the reaction rates; plasma physics 
effect enhance the radiative transition in the outgoing 
channel and enhance the reaction rates as well. 

The effects of neutron degeneracy on neutron capture 
reaction rates can be quantified by introducing the ratio 
R n , Eq. (U), of rates calculated for given conditions to 
those for non-degenerate neutrons. We have described 
this ratio by a simple analytic expression ([7]) assuming 
the power-law energy dependence of the reaction cross 
section §5§ at energies that are not too high. The de- 
rived expression ([7]) seems sufficient for many applica- 
tions. Furthermore, approximating cr(E) by a power-law 
function ([5]) , one also obtains the power-law index v and 
the maximum energy E max to which the power-law ap- 
proximation of cr{E) is valid. Eq and v are needed in 
Eq. ([7]), while E max controls the validity of Eq. ((Jj). 

Our conclusions are as follows: 

1. Neutron degeneracy can significantly affect (71,7) 
reactions in deep neutron star crust (Sec. HI]) . In 
many cases the effects of neutron degeneracy are 
well described by the factor R n given by Eq. ([JJ). 
For threshold reactions, strong neutron degeneracy 
enhances the reaction rate by many orders of mag- 
nitude. 

2. Plasma physics effects can additionally enhance 
(71,7) rates (Sec. IIII|) . that is described by the fac- 
tor R p \. These effects are less dramatic but can 
reach a few orders of magnitude. 

Furthermore, in Sec. II VI we have used the detailed bal- 
ance principle and calculated the rates of inverse (7,7*1) 
reactions taking into account neutron degeneracy and 
plasma effects. Finally, in Sec. [V] we discussed the ef- 
ficiency of (n.,7) reactions in a neutron star crust, with 
the conclusion that neutron degeneracy can be most im- 
portant. 

Finally it should be noted that free degenerate neu- 
trons in a neutron star crust can be in a superfluid state. 
Critical temperature T cn for the appearance of neutron 
superfluidity is very model dependent. Numerous cal- 
culations using different techniques (e.g., Ref. [3(J) give 
density-dependent T cn (p) with maximum values ranging 
from ~ 0.2 - 0.3 MeV to - 2 MeV, indicating that su- 
perfluidity is most likely. Superfluidity produces a gap 
in the energy spectrum of neutrons near the Fermi level 
and modifies matrix elements of neutron capture reac- 
tions. Both effects on neutron captures are not explored 
but may strongly modify the reaction rates. 

Our consideration of neutron degeneracy and plasma 
effects on neutron capture rates is simplified. A more 
rigorous (and complicated) analysis of these effects (in- 
cluding also neutron superfluidity) would be desirable. 
It would be instructive to perform self-consistent calcu- 
lations of the structure of atomic nuclei immersed in a 
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Fermi sea of free neutrons, taking into account a com- 
pression of the nuclei by free neutrons (e.g., Refs. [HE!)- 
For simplicity, we have used a model of free neutrons 
which occupy the space between atomic nuclei. At den- 
sities p not much higher than the neutron drip density, 
it is sufficiently accurate (as follows, for instance, from 
results of Ref. (HI). I R a self-consistent approach, this 
model should be replaced by a more elaborated unified 
treatment of neutrons bound in nuclei and free outside. 

In any case one should bear in mind that neutron cap- 
ture reactions in a deep neutron star crust can be af- 
fected by neutron degeneracy, plasma physics, and neu- 
tron superfluidity. These effects may have important im- 
pacts on nuclear burning and nucleosynthesis in the deep 
neutron star crust. The effects should be taken into ac- 
count to correctly simulate and interpret various obser- 
vational phenomena in accreting neutron stars such as 



X-ray bursts and supcrbursts as well as quiescent ther- 
mal emission of neutron stars in X-ray transients (e.g., 
Refs. [H, 0, H| and references therein). 
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Appendix: Approximation of Fermi-Dirac integrals 
by Aymerich-Humet et al. 

The asymptotes of Fermi-Dirac integrals (JSJ at y — > 
±oo in nondegenerate and strongly degenerate limits are 
easily obtained by retaining the first term in the series 
expansion of the Fermi-Dirac distribution function over 
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cxp(x - y): 



where the fit parameters are 



My) 
My) 



y -> -oo, 



v 



r(i/ + 2)- 



2/ 



(A.la) 
(A.lb) 



Aymerich-Humet et al. |27| derived a useful approxi- 
mation for Fermi-Dirac integrals, which, by construction, 
reproduces the asymptotic limits (|A.1[) . Their approxi- 
mation reads 



My) = 



r(^ + 2)2" +1 



[b + y + (\y - b\ c + a c )V l 



(A.2) 



a 



15, v 1 , 
l + -(v + l) + - Q (» 



b = 1.8 + 0.611/, 

c = 2 + (2 - v/2)2-". 



1/2 



(A.3a) 

(A.3b) 
(A.3c) 



In the range —0.9 < v < 4 for any y this approximation 
gives a relative error of « 1%. It is also valid for larger 
but become less accurate (to about 4% at v = 12). 



